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Intense, 100 fs laser pulses at 3.2 and 3.6 /im are used to generate, by multi-photon ionization, 
broadband wave packets with up to 400 eV of kinetic energy and charge states up to Xe^^. The 
multiple ionization pathways are well described by a white electron wave packet and field-free in- 
elastic cross sections, averaged over the intensity-dependent energy distribution for (e,ne) electron 
impact ionization. The analysis also suggests a contribution from a Ad core excitation in xenon. 
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Ionization in a strong laser field can be described by 
a three-step process IH [2] (Fig.l): production of photo- 
electrons, acceleration by the laser field, and (in)elastic 
recollisions or recombination with the parent ion result- 
ing in above-threshold ionization (ATI), high harmonic 
generation (HHG) or multiple ionization [3]-[6]. Each 
measureable conveys information on the intrinsic scat- 
tering event: a field-driven electron wave packet inter- 
acting with the core potential. Consequently, strong- 
field laser-driven scattering, or rescattering, is being vig- 
orously pursued as a method for imaging molecules [Tj. 
Recently, elastic rescattering was shown to yield informa- 
tion similar to conventional electron beam experiments 
[8^ , an established method for probing atomic and molec- 
ular structure [9 . The link between conventional and 
laser-driven scattering is expressed multiplicatively with 
a collisional cross section and a calculated electron return 
distribution that represents a quantum wave packet [10] . 
While both methods are believed to be comparable, laser- 
driven elastic scattering presents an advantage because it 
is synchronized to the period of the laser field. Conse- 
quently, it is considered a potential method for probing 
structure with sub- femtosecond resolution [TTl [12] . How- 
ever, as with all novel methods, laser-driven scattering 
requires more rigorous treatment that must be obtained 
through broader analysis. Currently, efforts for develop- 
ing laser-driven scattering have relied on the analysis of 
HHG and ATI. In this letter we investigate the third rec- 
ollision channel, inelastic scattering, by observing multi- 
ple ionization. Our experiment offers a new comprehen- 
sive basis for assessing laser-driven scattering and em- 
phasizes the use of inelastic events as a probe of complex 
electronic structure. 

Laser-driven inelastic scattering can result in non- 
sequential ionization (NSI) through (e,2e) [4l[5] or higher 
(e,ne) processes [6 , or excitation followed by field ioniza- 
tion [T3^. Most studies hitherto have been performed at 
near-infrared (NIR) wavelengths, (A < 1 fim). These in- 
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vestigations are limited to measurements of double NSI 
since depletion of the neutral ground state restricts the 
maximum intensity, and thus the return energy, to the 
(e,2e) inelastic channel. This limitation holds for all 
neutral atoms ionized by a NIR field. In addition, ex- 
traction of the inelastic NSI channel is masked by the 
presences of photon-assisted processes. In helium for ex- 
ample, semi-classical analyses [10 , 13 show contributions 
from both (e,2e) inelastic and field-ionized excited states. 
For higher- Z atoms, multiphoton [14] and higher-order 
sequential [151 US] processes complicate the analysis of 
the observed ion yields. Critically, the ability to exploit 
laser-driven scattering as a tool relies on the comprehen- 
sive understanding of all these experiments and the accu- 
racy of a field-free description of rescattering in a strong 
laser field. 

In this Letter, intense long wavelength pulses (3.6 and 
3.2 /im) result in a increase in the electron's rescat- 
tering energy [17] allowing the extension of NSI studies 
to (e,ne) processes for n=2-6. The corresponding obser- 
vations are both simplified and more comprehensive be- 
cause the field strength can be kept low enough to ionize 
only the neutral or excited singly charged ion and at the 
same time open a plethora of (e,ne) channels. We will 
show that our results can be interpreted by using known 
electron impact cross sections [18, 19 . Thus, the link be- 
tween NSI and laser-driven recollision is addressed over 
a broader, more significant, energy range. Moreover, at 
such impact energies, both inner-shell and valence ion- 
ization pathways are allowed in xenon [20 . This, by 
comparison, should illuminate the nature of laser-driven 
multiple ionization. 

We will show that laser-driven scattering is compa- 
rable to conventional electron impact ionization broad- 
ened by the nearly white returning wave packet that, 
according to classical mechanics, spans from to 3.2Up 
[TJ [2] where Up is the ponderomotive potential given by 
t/p[eV]=93Tntenstity[PW/cm2] • X^[j^m]. For compari- 
son, conventional experiments are performed with nearly 
monochromatic beams. A large Up and modest binding 
potential, /p, imply tunnel ionization since the Keldysh 
parameter [2P, 7 = y^Ip/2Up^ is less than unity. Our 
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FIG. 1: Diagram of the laser suppressed Coulomb potentials 
with the strong-field photoionization sequence illustrated: 
ionization, acceleration in the field with lateral expansion, 
and recollsion with a distorted potential. The Coulomb + 
laser potential of the neutral, dark blue solid line, and ion, 
light blue dashed line, are shown at a field intensity of 90 
TW/cm^ for the phase of birth that yields the maximum re- 
turn energy of 3.2Up. For reference the bare ionic Coulomb 
potential is shown by a thin black solid line. 



mid-infrared (MIR), 100 fs, laser is tuned to either 3.6 or 
3.2 /im and provides: 2b eV > Up >120 eV, 7 < 0.5, and 
electron wave packets with up to 400 eV of rescattering 
energy. 

The laser is described elsewhere ^2] and consists of 
a difference frequency generation (DFG) amplifier fed by 
two independent laser systems. The 100 fs pump beam is 
generated by a Ti: Sapphire chirped pulse amplifier and 
the 16 ps signal beam is generated by a regeneratively 
amplified NdiYLF laser [23^. DFG is performed in a KTA 
crystal and the idler beam has a central wavelength tun- 
able to 3.6 or 3.2 /im, a 100 fs pulse duration and a peak 
power exceeding 1 GW. The beam is focused into a time 
of flight mass spectrometer with a base pressure of 10~^ 
torr and 10^ — 10^^ cm~^ target densities. The mass 
resolution of the spectrometer is Am/m~0.3% and fully 
resolves all xenon isotopes up to the highest charge state 
observed. An intensity calibration, accurate to within 
±10%, was performed by measuring photoelectron spec- 
tra that have a characteristic change in slope ai 2Up [17]. 

Figure 2 shows the measured ion yields of xenon at 3.6 
jam where the highest charge state observed is Xe+^. At 
3.2 /im, not shown, the highest charge state observed is 
Xe+^. To verify that the ground state of the neutral is 
directly tied to the generation of higher charge states the 
sequential yield of Xe+^ was calculated using the Am- 
mosov, Delone and Krainov (ADK) [24 tunneling rate. 
Whereas the Xe+ measured yield is in close agreement 
with ADK for both 3.6 and 3.2 /im the predicted se- 
quential yield of Xe+^ underestimates the observed yield. 
In fact, Xe+^ is produced in greater abundance. Non- 
sequential behavior is well known in the NIR domain 
and should not be different in the MIR. However, at NIR 
wavelengths in xenon [Ml [15] non- sequential double ion- 
ization results from the multiphoton characteristics, e.g. 
7 > 1, of the ionization sequence and excitation. In the 
current study, observation of charge states up to Xe+^ is 
evidence of an extreme non-sequential process where over 
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FIG. 2: Ionization yields of xenon at 3.6 /im are shown with 
various symbols. The solid lines are calculated with ADK 
theory (see text). Top horizontal axis is the Keldysh param- 
eter. The vertical error bars are given by Poisson statistics 
and are, in general, smaller than the size of the symbols. Hor- 
izontal error bars show the 10% uncertainity in our intensity 
estimate. The calculated nonsequential yields are shown with 
dashed lines and, for reference, the ionization potentials for 
Xe^-Xe^^ are given. 



700 MIR photons are coupled to neutral xenon liberat- 
ing 6 electrons bound by a net potential of 243.6 eV. We 
have calculated non-sequential yields, shown with dashed 
lines in Fig. 2 and discussed below, with a nearly white 
wave packet and (e,ne) cross sections and found excellent 
agreement. 

Our analysis begins by defining inelastic branching ra- 
tios as Xe+"^/Xe+, n > 1, and assuming that the yield 
of the n^^ charge state is proportional to photoioniza- 
tion of the neutral. Since the returning electron flux is 
proportional to the ionization rate, our branching ratios, 
see Figs. 3(c) and 3(d), are effective, energy averaged 
cross sections for laser-driven multiple ionization. For 
comparison, the known experimental cross sections mea- 
sured with electron guns [18|il£^, are shown in Fig. 3(b). 
Experimental data is not, to the best of our knowledge, 
available for both quintuple ionization and 5s5p^ exci- 
tation of Xe+ so we have used the Lotz formula [25] to 
scale measured cross sections of Xe+. As an example, we 
observe a threshold response in our branching ratios that 
is similar to what occurs in the measured cross sections. 
This is most clearly seen in our experiment at 3.2 /im 
where the ponderomotive potential is 25% lower than at 
3.6 /im. 

As already mentioned, the spectral width of the return- 
ing electron distribution is known to be very broad, due 
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FIG. 3: The returning photoelectron distributions are shown 
for 100 (red), 80 (green), 64 (blue) and 51 (black) TW/cm^ in 
panel (a). For reference the coUisional excitation (5s5p^) and 
ionization cross sections used in our calculation are shown in 
panel (b), see text for details. The error bars for (e,2e)-(e,5e) 
are taken directly from refs. [181 119|. For 5s5p^ excitation 
and (e,6e) a Lotz formula is assumed with a minimum error of 
25%. Panels (b) and (c) show the results of Eqn. 1 for 3.6 and 
3.2 /im respectively, where the intensity axis is determined by 
the calculation, and the vertical overlap is referenced against 
the full NSI calculation shown in Fig. 2. 



to the distribution of ionization times {W . In addition, 
laser experiments are necessarily integrated over the spa- 
tial and temporal intensity distributions of the focussed 
beam. Hence, we define an effective, energy-averaged 
[26] , cross section, a(/) as: 
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where Wp is the net return distribution, and (j{E) is the 
experimental cross section. Consequently, the effective 
cross section is a function of the peak intensity /. 

To calculate Wp the laser pulse is sampled in time in- 
crements of 1 au and at each point a 1-D trajectory is 
calculated by solving x = —E{t); only trajectories that 
return to the core are retained. The initial conditions 
are consistent with the treatment in [IT, "W. Each tra- 
jectory is weighted by the ADK yield to account for tun- 
nel ionization and ground state depletion. The ballistic 
expansion of a tunnel ionized wave packet [27 is taken 
into account to ensure that each classical trajectory rep- 
resents a wave packet with the correct area at recollision. 
Finally, the calculation is summed for individual volume 
elements of a Gaussian focus. 

Examples of the return distributions are shown in Fig. 
3(a). Our classical representation of the returning wave 



packet strongly resembles the quantum mechanical ver- 
sion calculated in [10 . In general each distribution has a 
width of approximately 1.5/7p, increases slowly by a fac- 
tor of 3 to a maximum at around 2Up and then decreases 
rapidly at a maximum return energy of 3.2Up [1 . 

The average recolliding wave packet has an width 
of more than a nanometer. For reference, the (e,2e) 
cross section shown in Fig. 3(b) has a width of about 
1 Angstrom. Therefore, we need to include the depen- 
dence of the ionization probability on the impact param- 
eter [131 EH] since this determines the effective current 
density of the return distribution and quantifies the ob- 
served yield from the effective cross section Eqn. ([T]): 

exp(6Vag) 



P(6,/)=a(/)- 



(2) 



where ao = and b is the impact parameter. 

Our approach allows two analyses of laser-driven scat- 
tering: the first investigates the role of an energetically 
broad wave packet and relates a laser-driven measured 
branching ratio to a field-free measured electron impact 
cross section. The second accounts for the lateral expan- 
sion of the wave packet and quantifies the observed yield. 
For both cases we use the largest impact excitation cross 
section in xenon, 5s5p^, and (e,2e) impact ionization to 
account for Xe+^ yield, and (e,3e)-(e,6e) impact ioniza- 
tion cross sections for Xe+^-Xe+^, respectively. Xe+* has 
a binding potential of 9.9 eV and is assumed to ionize 
with unity probability. In Fig. 2, the yields, and in Figs. 
3(c) and 3(d), the branching ratios, from excitation and 
(e,2e) ionization are plotted separately to show the dif- 
ferent contributions from each channel. We see that for 
both analyses, and for both 3.6 and 3.2 fim laser fields, 
the yield predicted from excitation alone overestimates 
our experimental results while the observed multiple ion- 
ization channels agree with our assumption of a direct 
(e,ne) ionization event. We argue that excitation ion- 
ization in NSI may be beyond a semi-classical analysis 
because: first, excitation cross sections are typically in- 
ferred through a fitting procedure [28] and may have large 
errors. Second, recent work [29] suggests that higher ly- 
ing excited states are coupled to the field and may not 
be accurately described by field-free cross sections. 

Experimental and calculated NSI yields for Xe+^ and 
higher agree within a 10% uncertainty in our intensity es- 
timate and our overall experimental error (Fig. 2). In ad- 
dition, the excellent agreement between observed branch- 
ing ratios and effective cross sections for channels (e,ne) 
(n > 3) (Fig. 3 (c) and (d)) justifies the correspondence 
between laser-driven inelastic scattering and field-free im- 
pact ionization and emphasizes the key role of the ener- 
getically broad recolliding wave packet. Moreover, it jus- 
tifies retaining only channels that couple higher charge 
states to photoionization of the neutral. Finally, experi- 
ments performed with krypton for processes up to (e,5e) 
were found to be consistent with our interpretation. 

A qualitative argument that justifies our field- free in- 
terpretation appears in Fig. 1 which shows the Stark 
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FIG. 4: The measured (e,3e) cross section from [20^, red long 
dash, effective cross section from Eqn. 1, grey short dash, and 
the (e,3e) direct ionization channel, thin black line, obtained 
with a Lotz fit |25| are shown to highlight the feature between 
100-150 eV. Photoionization of a 4(i electron is indicated by 
the purple shaded region, see text for details. Measured +3/+ 
branching ratios for 3.2 /xm, green solid circle, and 3.6 /xm, 
green open circle, are shown on an energy scale determined 
by the peak of the return distribution that occurs at 2Up. 

ionic potentials for the maximum return energy of 3.2Up. 
Because the recolliding electron returns near the zero of 
the laser field the Coulomb potential suffers only a small 
distortion. We also note that even at the field maximum 
the region of the potential near /p(zon), light blue Gaus- 
sian in Fig. 1, is not significantly distorted. The success 
of the field-free cross sections has thus, a simple expla- 
nation in the three-step semi-classical model. 

Finally, the high ponderomotive energy available at 
long wavelengths opens more complicated ionization 



pathways. For example, (e,3e) ionization may involve 
an Auger (NOO) channel, or giant resonance, where 
the primary electron is removed from the M sub-shell 
[I9l |20]. Achenbach et al. [20 supported this argu- 
ment with a comparison of their electron impact ioniza- 
tion data with synchrotron photoionization experiments 
[30] . Their data along with the highlighted photoion- 
ization feature is shown in Fig. 4. For reference, both 
the effective, energy averaged cross section from Eqn. 1 
and the (e,3e) direct ionization cross section are shown 
in Fig. 4. Our (e,3e) data is consistent with the feature 
identified in [20 and has a small hump around 100 — 150 
eV that suggests inner-shell ionization is also occurring 
in the laser-driven experiment. The role of the energeti- 
cally broad wave packet here is that it averages out finer 
features. 

In conclusion we have shown that multiple ionization 
occurs from laser-driven scattering with a MIR laser. The 
maximum observed charge state is Xe+^ which is shown 
to occur through a direct field-free (e,6e) channel. The 
link between field-free impact ionization and laser-driven 
NSI is shown to be the broadband returning wave packet 
that effectively averages the cross sections. An inner-shell 
Auger process is suggested by the data and illustrates 
the effect of energy averaging. Our analysis of multiple 
ionization will play a critical role in helping to establish 
laser-driven scattering as a tool and an application of 
strong-field science. 
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